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In this paper we present a determinant quantum monte carlo study of the two dimensional Hub- 
bard model with random site disorder. We show that, as in the case of bond disorder, the system 
undergoes a transition from an Anderson insulating phase to a metallic phase as the on-site repul- 
sion U is increased beyond a critical value Uc. However, there appears to be no sharp signal of this 
metal-insulator transition in the screened site energies. We observe that, while the system remains 
metallic for interaction values up to twice Uc, the conductivity is maximal in the metallic phase just 
beyond Uc, and decreases for larger correlation. 
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Introduction 

The metal-insulator transition arising from the com- 
petition of randomness and interactions remains an in- 
triguing problem in condensed matter physics. For ex- 
ample, the question of the existence of a metallic phase in 
two dimensions, for which an experimental consensus had 
emerged in the lOSO'si has been revisited with new sam- 
ples over the last decade, with developments which have 
driven a considerable amount of new theoretical work?'^. 

Several interesting lines of study have emerged which 
explore the interplay of one-body potentials and two- 
body interactions in more general contexts. The 
superfluid-insulator transition has been studied in disor- 
dered, interacting boson systems, where the existence of 
a thermodynamic order parameter, the superfluid density 
Ps, as well as the greater ease of numerical simulations, 
has resulted in many definitive result a'^i^i^ . The coexis- 
tence of a metal and a Mott-Hubbard insulating phase 
in the disordered half-filled Hubbard model has been ex- 
plored using both numerical and analytical techniqueai. 
The existence of insulating phases away from commen- 
surate fermion filling has been explored in models with 
bimodal distributions of on-site chemical potential^. Fi- 
nally, the question of metallic phases arising from the ad- 
dition of correlations to a band insulator is drawing new 
attention^*'. 

The commonly cited qualitative picture of the ap- 
pearence of a metallic phase out of a disordered one is 
that the interactions act to screen the one-body poten- 
tial. While several Quantum Monte Carlo studies of dis- 
ordered interacting fermions exisliiiiiSii^ which demon- 
strate the possibility of a metallic phase, none have 
looked quantitatively at this screening in the Hubbard 
hamiltonian. 

In this paper, we will present results for the conductiv- 
ity and renormalized site energy of the two dimensional 
Anderson-Hubbard model, 

H = - {c]^ci„ + c\^Cja) 



+ U Xl"iT"ii + ^i'^j - A*)Kt + "ii) (1) 
j j 

Here Cj^,{cja) are fermion creation(destruction) opera- 
tors on site j for spin a and nja — Cja'^ja is the number 
operator, t is the hopping parameter, U the onsite re- 
pulsion, and /i and the global chemical potential and 
local site energies respectively. Each ej is drawn inde- 
pendently from a uniform distribution on (— iA,+iA). 
We choose t = 1 to set our scale of energy. 

Our key conclusion is that while increasing U can drive 
an Anderson insulating phase metallic, there appears to 
be no sharp signature of this transition in the variance 
of the renormalized site energies. This suggests that the 
metallic phase arises at least partially from an additional 
mechanism beyond a simple screening of the one-body 
potential. 



Numerical Approach 

We employ the determinant quantum monte carlo 
(DQMC) methodic. Since many descriptions of the ap- 
proach exist, we only provide a brief sketch here, focus- 
ing on those features most relevant to the present study. 
DQMC is an exact method to compute the properties 
of tight binding Hamiltonians on finite lattices. The 
inverse temperature /3 in the partition function is dis- 
cretized and an auxiliary ("Hubbard-Stratonovich") field 
introduced to decouple the interactions. The resulting 
quadratic form in fermion creation and destruction oper- 
ators is integrated out analytically, leaving a sum over the 
Hubbard-Stratonovich variables which can be performed 
stochastically. 

We have chosen the imaginary time discretization size 
small enough such that the systematic "Trotter" errors 
are comparable to the statistical errors associated with 
the monte carlo sampling and disorder averaging. Of 
greater concern in these simulations is the finite lattice 
size and, in particular, the possibility of a "false" signal 
of metallic behavior which would occur if the localization 
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length exceeds the lattice size. We have verified that in 
the phases we identify as metallic the localization length 
(computed at U ~ 0) is less than the lattice size. 

To investigate the metal-insulator transition, we look 
directly at the dc conductivity which we obtain from the 
current-current correlation function 



Me, r) = V (4,..c,, - 4c,+ . J ] 



(2) 



We compute the Fourier transform ja;(q, t) oijx{(-, t) and 
its correlation function KM^-^t) = 0:c(q, r) jx(-q, 0)). 
Using the general formalism of linear response theory, the 
dc conductivity is given by 



CTdc = lim^ 



ImAa.^(q = 0,Lj) 



(3) 



The frequency dependent conductivity is given by the 
integral transform 



duj 



TT 1 - e-'3'^ 



lmAxxiq,(^) (4) 



It is difficult to obtain ImAa:a;(q, w) by inverting this 
integral equation, because it requires the determination 
of Aa;a;(q, r) on a very fine mesh of imaginary times r 
with very high numerical accuracy. However, if we insert 
T = (3/2 in Eqn.^ the function multiplying ImA^^^; (q, t^) 
for low T effectively restricts the integral to small u, so 
that we may approximate ImAa;2.(q — 0, lu) by (TdcW- The 
frequency integral may now be evaluated analytically, 
leading to the result: 



A,,(q = 0,r = /3/2) 



(5) 



This approximation is expected to be valid when the 
temperature is smaller than an appropriate energy scale 
in the problem. It is convenient because it allows the 
computation of (Tdc as a function of temperature to be 
obtained from the function A^xi^l, t), which is calculated 
directly in DQMC. 

Obtaining a transport property like tJdc directly from 
imaginary time data, as described above, is a process 
which must be undertaken with caution. However, the 
use of this procedure gives the correct characterization as 
a metal or insulator in all cases which we have checked so 
far. For example, duiic/dT is positive (insulating behav- 
ior) for the half-filled d — 2 Hubbard model without ran- 
domness at all values of U , that is, regardless of whether 
the insulating character arises predominantly from an- 
tiferromagnetic order (weak coupling) or Mott behavior 
(strong coupling)iiii^. The procedure also gives the cor- 
rect physics in a band-insulator when a staggered site 
energy is present and J7 = 0. It has also been shown 
to give the correct physics of the disordered attractive 
Hubbard modeli^. 

A fundamental check of the numerical data is the veri- 
fication that the longitudinal current-current correlation 
function obeys the gauge invariance condition. 



where K, is the kinetic energy. We have checked that, as 
in previous workiiiiSiii this condition is satisfied. 

Since the disordered site-energies in the system 
are generated randomly from a uniform distribution 
(— ^A, -l-iA), the distribution has zero mean, and a vari- 
ance: 
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In order to study the screening of the disordered po- 
tential by interactions we note that within a mean-field 
picture an electron moving in the one body potential 
will feel the site energy renormalized by the density of 
oppositely oriented electrons. That is. 



UiUi 



which becomes, in the absence of spin-polarization, 

= (j + 



(8) 



(9) 
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since for each spin species {nja) — \{'nj). We define an 
associated dimensionless variance by normalizing to the 
fluctuations in the original site energies, 

= (^^)') - - (^~.)') (10) 

In the absence of interactions {U = 0), or for very large 
A at fixed U, we have V = 1, indicating that there is 
no screening of the random potential. The question we 
wish to address is whether there is some signal, e.g. a 
noticeable decrease, in V on entry into a metallic phase. 



Metallic phase due to interactions 

We begin by demonstrating that interactions can drive 
an Anderson insulating phase metallic. We show in Fig.Q] 
the dc conductivity as a function of temperature for a 
fixed strength of the disorder potential and increasing 
U. At U — 0, (Tdc decreases as T is lowered, indicating 
insulating behavior. However, at strong coupling (Jdc in- 
creases as T is lowered, indicating a cross-over to metallic 
behavior. All results in Fig. ^ and subsequently in this 
paper, are at one-quarter filling p = (n) — i. This is far 
away from the most dramatic effects of U in the Hub- 
bard model- the Mott transition and antiferromagnetic 
ordering. 

The metallic phase in Fig.^can be caused to return to 
insulating behavior by increasing the site disorder. This 
is shown in Fig. [21 where we begin with the interaction 
strength which gives the largest conductivity, U — A, 
and make A larger. For 9 < A < 10 the low temperature 
slope of fJdc reverts to insulating character. 

An interesting feature of Fig. ^ is the non-monotonic 
behavior of the conductivity, ddc increases with U up to 
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FIG. 1: (Color online)Tlie dc conductivity as a function 
of temperature for increasing values of the on-site repulsion 
U = — 5. The site energy variance A = 8. 
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FIG. 3: (Color online)Like Fig.Q except at larger disorder, 
A = 9. The same decrease of conductivity with U in the 
metallic phase is seen as in Fig.0 
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FIG. 2: (Color online)The dc conductivity as a function of 
temperature for increasing values of disorder A — 8 — 11. The 
on-site repulsion 17 = 4. 



U ~ 3 — A, but then comes down again at U — 5. In or- 
der to verify that this phenomenon is generic, we show in 
Fig.|3|data for larger A = 9. We again see that a^c comes 
down at strong coupling. A similar phenomenon occurs 
in the evolution of the superfluid density ps for corre- 
lated bosons moving in a random potential- a superfluid 
phase with ps ^ exists at intermediate coupling, but 
the system is insulating, ps — 0, both at weak and strong 
coupling^ . 

In Fig. 2] we show a^c vs. the disorder strength for 
progressively lower temperature values. For A < Ac, the 
system is metallic and adc increases as the temperature 
is lowered, while for A > Ac, in the insulating state, 
the behavior is opposite. The crossing point of the plots 
demarcates the critical disorder strength. 

In Fig. [SI we show a similar crossing plot for ddc as we 
tune the interaction strength through the metal-insulator 



FIG. 4: (Color online)A crossing plot for a^c vs. A. The 

critical disorder strength Ac ~ 9.2 — 9.3 for (/ = 4 is clearly 
seen. 
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FIG. 5: (Color online)A crossing plot for (Jdc vs. U. 
crossing is seen to happen for < U < 1 for A = 8. 
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transition for a fixed disorder strength A = 8 (See Fig.^. 
A small value of the interaction is seen to be enough to 
cause the transition to a metal. Interestingly, the con- 
ductivity is non-monotonic and decreases for large values 
of the interaction strength (The fermion sign problem 
in DQMC simulations forbids the evaluation of adc at 
U — = 8). It is possible that there is a crossing 
at larger interaction strengths, when the system reverts 
back to an insulator. Such non-monotonic behavior be- 
havior of the conductivity has also been seen in recent 
DQMC studies of a multi-band Hubbard model at half- 
filling , where the sequence of transitions with increasing 
U is found to be Band insulator metal Mott insu- 
lator. 



Renormalized site-energies with interactions 

Before showing results for the dimensionless variance of 
the renormalized site energies, we present in Fig.|Sla plot 
of the original and renormalized site energy landscapes. 
As expected, the particles preferentially sit on the sites 
with low Ej , and these larger values of {rij ) then lead to a 
smoother ej = Cj + ^{nj). However, there is certainly no 
very dramatic leveling of the landscape. Below, we will 
explore this more quantitatively. 
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FIG. 6: (Color online)Left: Landscape of the original site 
energies tj with A = 8. Right: Landscape of the renormalized 
site energies ij with A = 8, U = 4 and /3 = 8. On the right, 
the mean increase in the renormalized site energies due to U 
has been subtracted out. 

In Fig. \7\ we examine whether there is signal of the 
metal insulator transition in the evolution of V. We 
plot the low temperature slope dadc/dT from the data 
of Fig. 121 and show its change of sign at A « 9.2 — 9.3. 
There is no clear indication of this critical value in the 
renormalized site energy variance V. 

We can similarly look for this effect at the metal- 
insulator transition driven by increasing U at fixed A = 8 
(Fig. This is shown in Fig. |S1 Again, there appears 
to be no clear signal of the metal-insulator transition in 
the screened site energies. 




FIG. 7: (Color online)The variance of the renormalized site 
energy V is shown as a function of A, as is the low temperature 
slope of the conductivity {U = 4, (3 — 8) . 
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FIG. 8: (Color online)The renormalized site energy is shown 
as a function of U. There appears to be no signal of the MIT 
at small U, nor the conductivity peak at ?7 ~ 3 — 4. 



Renormalized site-energies with zero hopping 

The results from the previous section suggest we look 
more closely at the physical picture of the smoothening 
of the site-energy landscape by interactions. Our expec- 
tation in Fig. Ul where we plotted V as a function of the 
disorder strength was that in the metallic phase at weak 
disorder there would be a markedly smaller value of V, 
and then a crossover to a larger value as the disorder is 
increased into an insulating phase. On the other hand, 
at weak disorder we expect the least inhomogeneity in 
the site occupations. In the limit of uniform density, at 
very weak disorder, the site energy variance equals the 
original one, and V'^ — 1, suggesting that V might in- 
stead decrease with disorder. The conflicting tendency 
to decrease with disorder as charge-inhomogeneity devel- 
ops, and increase in the insulating phase, might explain 



5 



why the site energy variance is so insensitive to site en- 
ergy disorder, whereas when we tune through the transi- 
tion with interaction strength there is a much larger de- 
chne (though still no abrupt signal at the transition). In 
this section, we examine an analytically solvable limit of 
the disordered Anderson-Hubbard model, that of t = 0, 
which can be considered to be the limit of very high dis- 
order. At this limit, there is no metallic behavior, but 
it is still interesting to investigate the behavior of the 
site-energy distribution as we move from weak to strong 
interaction (or, strong to weak disorder). 
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FIG. 9: (Color online)At t = the energy levels for occu- 
pation by the first electrons (which we denote by j) extend 
from —A/2 to A/2 (left). The energy levels for occupation 
by the second electrons (which we denote by J,) extends from 
-A/2 + U to +A/2 + U (right). 



When there is no hopping, the Anderson-Hubbard 
model is classical. As electrons are added to the lat- 
tice at t = 0, the sites with the lowest site energies are 
singly occupied up the Fermi energy Ep. When, how- 
ever, Ep exceeds —A/2 + U it becomes preferable to 
start doubly occupying the low energy sites. This is il- 
lustrated in Fig. 1^1 From the figure it is evident that 
(n^) = ^{Ep + A/2) and (n^) = ^{Ep + A/2- U), and 
hence that (n) = {n-f + ni) = ^{2Ep + A - [/). We 
can easily obtain the mean of the renormalized site ener- 
gies by averaging Eqn.|Hl (e|) + (q) = U{(n^} + (nj,)) — 
U{2Ep + A-U). 



A completely equivalent result is obtained by recog- 
nizing that the energies of the sites on which up spin 
electrons reside are raised by U in the range from — A/ 2 
to Ep — U, where down spins are present. Similarly, the 
energies of the sites on which down spin electrons reside 
are raised by U in the range from — A/2 to £'f, where up 
spins are present. (Here the designations 'up' and 'down' 
merely reflect the 'first' and 'second' electrons on a site.) 
When the energies are averaged over these ranges, the 
same result for (e|) -|- (q) is obtained. 

The average of the square of the renormalized site en- 
ergies is obtained in the same way. We can then evaluate 
the dimensionless variance of the renormalized site ener- 



defined as 



12 
A2 



(11) 



To determine the variance, we must distinguish between 
two cases: a generic one, in which the Fermi energy Ep 
is larger than — A/2 -I- [/, and there is double occupancy 
of the low-energy sites, and a non-generic case, in which 
there is only single occupancy, which may happen for a 
large value of x = J7/ A or a small density. For the generic 
case, the variance can be computed in terms of the three 
energy-scales (J7,A and Ep) in the t = problem: 
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(12) 



The Fermi energy (for the generic case) can be deter- 
mined in terms of U and A by the following equation: 



P = 



2 ^ 



(13) 



where p is the filling. For example, in the quarter-filled 
case, p = i, and Ep = ^{U ~ y)- l'^ the non-generic 
case, the Fermi energy can, of course, be determined by 
simple state counting. 

A plot of Vf^Q vs. X = U/A is given in Fig. ^ for dif- 
ferent fillings. At very weak interaction (or, very strong 
disorder), the variance equals the non-interacting value 1. 
As the interaction is increased (equivalently, the disorder 
is decreased), the variance first decreases and reaches a 
minimum. Upon further increasing the interaction, how- 
ever, the site occupations become homogeneous and the 
variance grows. 



Conclusions 

We have examined the metal-insulator transition in the 
Anderson-Hubbard model using determinant Quantum 
Monte Carlo simulations. Our focus has been on the evo- 
lution of the renormalized site energy through the tran- 
sition, and we conclude that it exhibits no sharp feature 
there. It seems that the picture of screening of the disor- 
der by interaction is too primitive to account for metallic 
behavior. 

On the other hand, we observe an interesting non- 
monotonic behavior of the conductivity with interaction 
strength. In the boson Anderson-Hubbard model, the 
ground state at incommensurate densities is an Ander- 
son insulator at weak U and a insulating "Bose glass" at 
large U. In between, there is a superfluid phase in which 
the superfluid density first rises as one emerges from the 
Anderson insulator and then falls to zero again upon en- 
try into the Bose glass. As far as we can see, the fermion 
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Hubbard model remains metallic at large C/, but simula- 
tions there are difficult and we cannot make a definitive 
statement. In any case, the non-monotonic behavior of 
(Tdc is rather analogous to the behavior seen for strongly 
interacting, disordered Bose systems. 



FIG. 10; The site-energy variance, as defined in Eq. 1111 is 
plotted for three different values of the filling, p — 1, 0.5 and 
0.25. The x-axis is the ratio of the interaction to the disorder 
strength, U/A. 
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